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Abstract. Over the past six years, the Fermi Large Area Telescope has detected more than 
150 7 -ray pulsars, discovering a variety of light curve trends and classes. Such diversity hints 
at distinct underlying magnetospheric and/or emission geometries. We implemented an offset- 
dipole magnetic field, with an offset characterised by parameters epsilon and magnetic azimuthal 
angle, in an existing geometric pulsar modelling code which already includes static and retarded 
vacuum dipole fields. We use these different R-field solutions in conjunction with standard 
emission geometries, namely the two-pole caustic and outer gap models (the latter only for 
non-offset dipoles), and construct intensity maps and light curves for several pulsar parameters. 
We compare our model light curves to the Vela data from the second pulsar catalogue of Fermi. 
We use a refined chi-square grid search method for finding best-fit light curves for each of the 
different models. Our best fit is for the retarded vacuum dipole field and the outer gap model. 


1. Introduction 

The discovery of the first pulsar in 1967 by Bell and Hewish m gave birth to pulsar astronomy. 
Pulsars are believed to be rapidly-rotating, compact neutron stars that possess strong magnetic, 
electric, and gravitational fields [T]. They emit radiation across the entire electromagnetic 
spectrum, including radio, optical, X-ray, and y-rays [ 1 ]. Since the launch of the Fermi Gamma- 
ray Space Telescope in June 2008, over 150 y-ray pulsars have been detected, of which the Crab 
and Vela pulsars are the brightest sources. Fermi consists of two parts including the Large Area 
Telescope (LAT) and the Gamma-ray Burst Monitor. The LAT measures 7-rays in the energy 
range between 20 MeV and 300 GeV [3j. Over the past six years, Fermi has released two pulsar 
catalogues, both describing the light curve profiles and spectral characteristics of 7-ray pulsars 
m- The light curves show great variety in profile shape, and may be divided into three general 
classes based on the relative phase differences between the radio and 7-ray pulses [IHIM]. The 
light curves also show energy-dependent behaviour. Most of the young and millisecond pulsars 
exhibit two 7-ray peaks whereas some pulsars including Vela display three peaks [ 2 ]. 

Models of pulsar geometry are characterised by the inclination angle (a) between the rotation 
(fi) and magnetic (fx) axes, the observer’s viewing angle (C) between the rotation axis and the 
observer’s line of sight, and the impact angle (/3 =| </ — a |). Geometric models assume the 



Figure 1. A schematic representation of the 
different geometric pulsar models. The PC 
model extends from the neutron star surface 
^ up to low-altitudes above the surface (yellow 
region). The TPC emission region (curved 
magenta lines) extends from R^s (neutron star 
radius) up to the light cylinder i?LC (where 
corotation speed equals the speed of light, 
vertical lines), the OG region (cyan regions) 
from i?NCS (the null charge surface, where the 
Goldreich-Julian charge density pcj = 0 m, 
indicated by the dark blue lines) up to the RhC, 
and the PSPC from i?NS to the Rlc, covering 
the full open volume region. Adapted from [El- 


presence of several ‘gap regions’ in the pulsar magnetosphere. These are defined as regions 
where particle acceleration and emission take place. These geometric models include the two- 
pole caustic (TPC) |9] (the slot gap (SG) [19] model may be its physical representation), outer 
gap (OG) [ 6 l | 22 ], and pair-starved polar cap (PSPC) model [15]. All of these models are 
represented in hgure [TJ The emissivity Cj/ of high-energy photons within this gap region is 
assumed to be uniform in the corotating frame for geometric models (for physical models the Cj/ 
changes with radial distance when an electric field is assumed). Since the 7 -rays are expected 
to be emitted tangentially to the local magnetic field in the corotating frame m, the assumed 
magnetic held geometry is very important with respect to the predicted light curves. Several 
magnetospheric structures have been studied, including the static dipole held m, the retarded 
vacuum dipole held (RVD) [ 8 ] and the offset-dipole R-held. The latter is motivated by the fact 
that retardation of the R-held at the light cylinder causes offset of the polar caps (PCs), always 
toward the trailing edge of the PC {4>o = 7 r/ 2 ) [16] . The offset is characterised by parameters 
epsilon (e) and magnetic azimuthal angle (0o) which represents a shift of the PC away from the 
magnetic axis, with e = 0 corresponding to the static-dipole case. 

In this paper we studied the effect of using different combinations of magnetospheric 
structures, geometric models, and model parameters on 7 -ray light curves. In Section [2] we 
discuss the implementation of an offset-dipole solution and associated electric held. Section [3| 
describes how we matched limiting cases of the low-altitude and high-altitude F-helds using 
a matching parameter rjc- In Section H) we describe the chi-squared (y^) method we used to 
search for best-ht light curves. Our results are given in Section [5] and the conclusions follow in 
Section [ 6 ] 

2. Implementation of an offset-dipole R-field 

We implemented [5] an offset-dipole magnetic held m into an existing geometric pulsar 
modelling code [10] which already includes static and RVD helds. The implementation involves 
a transformation of the R-held from spherical to Cartesian coordinates, rotating both the B- 
held components and its Cartesian frame through an angle —a, thereby transforming the R-held 
from the magnetic frame (z' || p) to the rotational frame (z || ft). We extended the range of 
e for which we could solve the PC rim (for details, see [TO]) by enlarging the colatitude range 
thought to contain the last open held line (tangent to Rlc, see hgure [T]). 















3. Matching parameter 

It is important to take the accelerating ill-field into account (in a physical model) when such 
expressions are available, since this will modulate the in the gap (as opposed to geometric 
models where we just assume constant at all altitudes). We use analytic expressions [20] for 
a low-altitude and high-altitude SG £'-held in an offset-dipole magnetic geometry. These are 
matched at a critical scaled radius rjc = rd i?NS to obtain a general ill-field valid for all altitudes 
|2Uj . In previous work we chose r]c = 1.4. In this paper we solve r]c on each field line. Using 
the general ill-field we could solve the particle transport equation [231E] (taking only curvature 
losses into account) to obtain the particle Lorentz factor 7 that is necessary for calculating the 
curvature emissivity. 


4. Finding best-fit light curves 

We tested several fitting methods in order to select the most suitable one. We decided to use 
the standard method. For each combination of ii-field and geometric model, we calculate 
for each set of free model parameters a, C, normalisation (A) and phase shift (A(/>). We 
assume that for the bright Vela pulsar, with a large amount of counts in each bin, the x^ will 
be Gaussian distributed with N^of = 96 the number of degrees of freedom. For the standard 
Gaussian distribution, assuming very small Gaussian errors, the reduced Xu values are very large 
and therefore we needed to scale the x^ values with the optimal Xopt ™cl multiply by N^oi- The 
scaled x^ is denoted by and the reduced scaled x^ will have a value of e = im- We next 
determined confidence intervals (Icr, 2a, and 3a) in (a,C) space for these model parameters, 
using Vdof = 2. We only fitted the 7 -ray light curve, because we do not want to bias our results 
with a simplistic radio model. 

5. Results 

5.1. Phaseplots and light curves 
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Figure 2. Phaseplots and light curves for Figure 3. Phaseplots and light curves for 
different e values, for constant Ci,. different e values, for variable 



























C =90° 



Phase iti 




1 \ 

1 \ 


/ 


,1 1, 

f \ 

i> \ 

/ 

f \ 

J \ 

A 

/ \ 



1 r 


f\ 

/ \ 

/ \ 

i', 

11 

vv 


A , 
, “ \ ' 
i / 


; V 

'1 

l\ 

' \ ./V-, 

\ /\ / 
V '■'*7 

K 


f\hj 

1: 1 V 

) A, 

1 

f- 

n. 

jl 

1 ^ 


1 1 

i i 

1 

1 

1 ,1 




i 

. , . .1 L 

1 



0 • • • • .... .... ... .1 Li-----■- . .—^-.-1 

0 Q£ OA 0.6 D.S 0.2 0.4 0.6 0.8 0.2 04 0.6 0.8 0.2 04 0.6 O.S 0.2 0.4 0£ 0.8 0.2 0.4 0.6 0.8 1 

Phase 4* 


Figure 4. Phaseplots and light curves for different values of a and Ccut for a hxed value of e = 0.18, 
for variable Cj/. 


In figure [21 [3] and m we show the intensity maps or phaseplots (emission per solid angle 
versus ( and pulse phase 4>) and their corresponding light curves (i.e., cuts at constant C = Ccut) 
for the offset-dipole S-field and the TPC model. The dark circle in figure [2] is the non-emitting 
PC, and the sharp, bright regions are the emission caustics, where radiation is bunched in phase 
due to relativistic effects. Figure[2]and [3]represent phaseplots for a fixed a = 40° and Ccut = 70°, 
with e ranging from 0 to 0.18 with increments of 0.03. We contrasted the cases of constant and 
variable Ey. We observed a qualitative difference in caustic structure. The caustics seem larger 
and more pronounced in the constant case. In figure [4] we chose a fixed value of e = 0.18 for 
variable e^/, with a ranging from 0 ° to 90° and Ccut from 15° to 90°, both with a resolution of 
15°. This shows examples of various light curves that may be obtained in this model. 

5.2. Contours and best-fit light curves 

In figure [5] we represent our best fits we obtained using the method. Our overall best fit 
is for the RVD field and the OG model, with a = 78°, C = 69°, A = 1.3, and Ac/) = 0. For 
the offset-dipole solution and the TPC model, we have a best fit (assuming constant e^) for 
parameters e = 0.00, a = 73°, C = 45°, A = 1.3, and Ac/) = 0.55. The best fit for the offset- 
dipole solution and the TPC model, assuming variable Cy, is for parameters e = 0.18, a = 73°, 
C = 17°, A = 0.5, and Ac/) = 0.60. The best-fit parameters for each R-field and geometric model 
combination are summarised in table [TJ The table includes the different model combinations, 
the optimal yfi value (before scaling), the free parameters with errors (found using 3 cj connected 
(ca,C) contours), a reference fit found using radio polarisation data, and the comparison between 
models using the difference between the respective optimal values of represented by AC*. 

6. Conclusions 

We have studied the effect of implementing the offset-dipole R-field on 7 -ray light curves for 
the TPC geometry. We find an optimal best-fit for Vela for the RVD R-field and the OG 
model. The OG model displays reduced off-peak emission. We note that the best fits for the 
offset-dipole R-field for constant Cy favour smaller values of e and for variable Cy larger e values. 
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Figure 5. Contour plots (left) and corresponding best-fit light curves (right). Panels (a) and 
(b) represent the best ht for the RVD and the OG model. Panels (c) and (d) are for the offset 
dipole and TPC model for constant ei, and e = 0.00. Panels (e) and (f) are for the offset dipole 
and TPC model for variable €iy and e = 0.18. The colour bar of the contour plots represent 
logioC^. The conhdence contour for la (magenta line), 2cr (green line), and 3a (red line) are 
also shown. The star indicates the best-ht solution. The blue histogram denotes the observed 
Vela prohle (for energies > 100 MeV) [2] and the red line the geometric model. 


When including an F-held the resulting phaseplots becomes qualitatively different compared to 
constant €i,. In future, we want to continue to produce light curves using improved geometric 
models and R-helds, and also using more data, in order to search for best-ht prohles, thereby 
constraining the low-altitude magnetic structure and system geometry of several bright pulsars. 
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Table 1. Best-fit parameters for each B-field and geometric model combination. 
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